clear;clc
line1='AGTACAAATAGCCGACGGCATTAA'%初始DNA链
A=strfind(line1,'A');%查找A字符的位置
T=strfind(line1,'T');%查找T字符的位置
C=strfind(line1,'C');%查找C字符的位置
G=strfind(line1,'G');%查找G字符的位置
n=length(line1);%统计初始链的长度
SUM=1:n;%生成第二条链
SUM(A)='T';
SUM(T)='A';
SUM(C)='G';
SUM(G)='C';%对第二条链的相应位置赋值对应碱基
line=char(SUM)%将 数字转化为字符
num=length(strfind(line1,'T'))%T-A碱基对的个数
